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ABSTRACT 

We present first results from the Palomar Adaptive Optics Survey of Young Stars conducted at the 
Hale 5 m telescope. Through direct imaging we have discovered a brown dwarf and two low-mass 
stellar companions to the young solar-type stars HD 49197, HD 129333 (EK Dra), and V522 Per, and 
confirmed a previously suspected companion to RX J0329. 1-1-0118 (Stcrzik et al. 1997), at respective 
separations of 0.95" (43 AU), 0.74" (25 AU), 2.09" (400 AU), and 3.78" (380 AU). Physical association 
of each binary system is established through common proper motion and/or low-resolution infrared 
spectroscopy. Based on the companion spectral types, we estimate their masses at 0.06, 0.20, 0.13, 
and 0.20 Mq, respectively. From analysis of our imaging data combined with archival radial velocity 
data, we find that the spatially resolved companion to H P 129333 is potentia.llv id entical to the 
previously identified spectroscopic companion to this star ijDuauennov fc Mavoilll99l|l . However, a 
discrepancy with the absolute magnitude suggests that the two companions could also be distinct, with 
the resolved one being the outermost component of a triple system. The brown dwarf HD 49197B is 
a new member of a growing list of directly imaged sub-stellar companions at 10-1000 AU separations 
from main sequence stars, indicating that such brown dwarfs may be more common than initially 
speculated. 

Subject headings: binaries: close — instrumentation: adaptive optics — stars: individual (HD 49197, 
HD 129333, V522 Per, RX J0329. 1+0118) — stars: low-mass, brown dwarfs 



1. INTRODUCTION 

High-contrast imaging searches for low-mass compan- 
ions to nearby and/or young stars have increased dra- 
matically in number since the initial discovery of a 
brown dwarf companion to a main sequence st ar (Gl 229) 
through direct imaging (Naka nma et a] .11 9951) . One par- 
ticularly powerful technique is adaptive optics (AO), 
which provides the high angular resolution (< 0.1") 
achievable at the diffraction limit of large ground-based 
telescopes. The widening niche of high-contrast imaging 
opened by recent developments in AO technology implies 
that not only brown dwarfs but exo-solar planets may be 
within the realm of direct imaging. Nowadays nearly ev- 
ery ground-based telescope equipped with an AO system 
hosts an imaging companion-search project. The sud- 
den explosion in interest in this topic has been fueled 
by the success of the radial velocit y (r.v.) method in de - 
tecting solar system analogs fe.g.. lMarcv fc Butlerll99'^ . 
Through longer time- lines of observation and higher pre- 
cision, the sensitivity of r.v. surveys has now extended 
outwards to include planets at semi-major axes >3 AU 
(jCarter et al. 2003: Naef ct al. 2004), i.e., near the Jo- 
vian region in the Solar System. While the sensitiv- 
ity of direct i maging to "planetary -mass" (1-13 Jupiter 
masses (^Mr): iBurrows et al.lll997|) objects at such sep- 
arations from Sun-like stars is still extremely limited 
due to contrast requirements, several higher-mass brown 
dwarf companions have been discovered at wider separa- 
tions (see compilation in lReid et al.ll2001tl. some at pro- 
jected distance s as small as 14-19 AU l)Els et al.ll200lt 
iLiu et al1l2002t) . High-resolution spectroscopic monitor- 
Electronic address: [metchev, lah@astro.caltech.ediil 



ing and direct imaging are thus complementary in search- 
ing for sub-stellar companions to stars. Future develop- 
ment of both methods promises to narrow, and eventu- 
ally close the sensitivity gap between them. 

Young stars are the most suitable targets for direct 
imaging of sub-stellar companions. At ages of 10- 
100 million years (Myr) the expected brightness ra- 
tio in the near IR between a 10 Mjup objec t and 
a so lar- type star is 10~^-10~^ ijBurrows et al.l Il997t 
Bara ffe et alJ I2003D . Modern AO systems can rou- 
tinely achieve comparable dynamic range at 1" sep- 
arations from bright stars. Hence, for young stars 
within 40 pc of the Sun we can probe for massive plan- 
ets at separations comparable to the giant planet re- 
gion in our own Solar System. However, few young 
stars are known at such small heliocentric distances. 
These are constrained to several tens of members of 

fung moving groups: TW Hya l|R,ucinski fc Krautteil 
83; dc La Reza et al. 1989), Tucana/Ho r ologium 
uc kerman fc Webbn'2 000: Zuckcrman et a ll [2001atl . 
and d Pic l|Zuckerman. Song, fc Webb .2001bE and 
have already been targeted with sensitive space-based 
and ground-based AO surveys, which have uncov- 
ered 3-4 brow n dw arf companions (Lowrance 2002 
ILowra,nce et 311^2003': -Neuhauser fc CuentheHl2004f but 
no planetary-mass ones. 

Because contrast and projected separation are the 
limiting factors in detectability of sub-stellar compan- 
ions, brown dwarfs, being more luminous than planets, 
are detectable at greater heliocentric distances and at 
smaller angular separations from their host stars. At 
the same physical separation from the primary (e.g., 
50-100 AU), brown dwarf companions should be de- 
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tectable around older (several gigayears [Gyr]) and/or 
more distant (^200 pc) stars compared to planets, al- 
lowing a larger sample of targets. With regard to this, 
we have commenced a survey of young (<400 Myr) 
solar-type (F5-K5) stars within 160 pc using the AO 
system on the Palomar 5-m telescope. Our survey 
sample is largely a subset of the sample targeted by 
the Formation and Evolution of Planeta ry Sy stems 
(PEPS) Spitzer Legacy Team ijMever et al.ll2004 . Al- 
though faint primary stars, such as M dwarfs or white 
dwarfs, offer more favorable contrast for imaging sub- 
stellar companions, we have chosen to limit our sam- 
ple to solar analogs because of interest in determining 
the multiplicity statistics of sub-stellar objects around 
other suns. Furthermore, several recent large sur- 
veys have already exp lored the multipli c ity of nearby 
(<50 pc) cool stars llClose et al.l 120031 iCarson et al.l 
I2003t iMcCarthy fc Zuckerman|l200 4D, or white dwarfs 
JZuckerman & Becklin 1992; Gree n. Ali. &: Napiwotzka 
12000), while a large sample of F-G stars has not been 
studied, because of comparatively small numbers in the 
immediate solar neighborhood. 

Preliminary results from our su r vey w ere reported 
in iMetchev. Hillenbrand, fc MeveJ lj2002(l . Here we 
present the strategy of the survey, and the discov- 
ery and confirmation of resolved low-mass companions 
to HD 49197, HD 129333 (EK Dra), V522 Per, and 
RX J0329. 1-^0118. We shaU refer to these through- 
out the paper as HD 49197B, HD 129333B, V522 PerB, 
and RX J0329.1-I-0118B. For convenience of notation, 
a second candidate companion to HD 49197 found to 
be an unrelated background star will be denoted as 
HD 49197"C". The full sample and further results from 
the survey will be discussed in a later paper. 

2. OBSERVING STRATEGY 

The observations described in this Section are repre- 
sentative of our general survey observing strategy. Ta- 
ble ^ details the imaging and spectroscopic observations 
specifically for the four objects presented here. The prop- 
erties of the observed primaries are given in Table |21 

2.1. Imaging 

2.1.1. First-Epoch Imaging and Survey Sample 
Subdivision 

First epoch observations a re obtained with the Palo- 
mar AO system (PAL AO; iTrov et all 12000) in resi- 
dence at the Cassegrain focus of the Palomar 5-m tele- 
scope. Since the summer of 2003 the wavefront sen- 
sor runs at frame rates up to 2 kHz, and the system 
routinely produces diffraction-limited images (0.09" at 
Ks) with Strehl ratios in the 30-50% range at 2/xm on 
V < 12 guide stars, and up to 75% on F < 7 mag 
guide stars. PAL AO employs the Palomar High Angular 
Resolution Observer (PHARO; Havward et al. 2001), a 
1024x1024 pix HgCdTc HAWAH detector with imaging 
(25 mas/pix and 40 mas/pix plate scale) and spectro- 
scopic (i?=500-2000) capabilities in the near IR. A set 
of coronagraphic spots, Lyot masks and neutral density 
(ND) filters are available to achieve the desired dynamic 
range. 

Our program entails ifs-band (2.15/im) imaging in the 
25 mas/pix mode (25" x 25" field of view) both with and 



without a 0.97"-diameter coronagraphic stop. For high 
dynamic range, long (1 min) coronagraphic images are 
taken to identify fainter (potentially sub-stellar) compan- 
ions at separations >0.5". Twenty-four such exposures 
are taken, for a total of 24 min integration per target, 
with 6 min spent at each of 4 different orthogonal de- 
tector orientations (obtained by rotating the Cassegrain 
ring of the telescope). For every 6 min of on-target imag- 
ing (i.e., at each detector orientation), separate 1-min 
coronagraphic exposures are taken at five dithered sky 
positions 32-60" from the star. For high angular resolu- 
tion (but with lower dynamic range), short (1.4-9.8 sec) 
non-coronagraphic exposures are taken to look for close 
companions of modest flux ratio, and to establish rela- 
tive photometric calibration. The images are taken in 
a 5-point dither pattern at the vertices and center of a 
box 6" on a side. A 1%-transmission ND filter is used if 
necessary to avoid saturation^. On occasion, a narrow- 
band (1%) Brackett-7 (2.17/im) filter is used for higher 
throughput, instead of the ND 1% filter. 

To avoid detector saturation and/or decreased sensi- 
tivity over a substantial fraction of the image area, stars 
with bright (Aifg < 4) projected companions in the 
PHARO field of view (FOV) were not observed with deep 
coronagraphic exposures. However, binaries with sepa- 
ration < 0.5" were included, as both components of the 
binary could then be occulted by the coronagraph. This 
naturally splits our survey sample in two groups: the 
"deep" subsample, consisting of essentially single stars 
and close binaries, and the "shallow" subsample encom- 
passing the remaining stars. Membership to one of the 
two subsamples was assigned at the telescope, when their 
multiplicity and approximate flux ratio was revealed dur- 
ing the short exposures. The shallow subsample was fur- 
ther expanded to include stars out to 200 pc and/or older 
than 400 Myr to cover the entire FEPS sample accessible 
from the Northern hemisphere. 

Short dithered exposures were taken of all stars, while 
long coronagraphic exposures were taken only of stars in 
the deep subsample at Ks band. In addition, short J- 
(1.22/im) and _ff-band (1.65/im) exposures were taken of 
all candidate binaries (all stars in the shallow survey, and 
the < 0.5" systems in the deep survey) to allow approx- 
imate photospheric characterization of the components. 

In accordance with the above distinction, HD 49197 
was observed for a total of 24 min with the corona- 
graph as a part of the deep survey, while HD 129333 
and RX J0329. 1-1-0118 (with bright candidate compan- 
ions), an d V522 Per (a P er member, 190 pc from 
the Sun; Ivan LeeuwenI 1199'^ were observed only with 
short exposures. Conditions were photometric during 
the first epoch observations of HD 49197, V522 Per, 
and RX J0329. 1-1-0118, and unstable during those of 
HD 129333. 

2.1.2. Follow- Up Imaging 

We obtain second-epoch imaging observations of all 
candidate companions to check for common proper mo- 
tion with their corresponding stars. Such were taken for 

^ The ND 1% filter was caUbrated photometrically through 
repeated (17-20 per band) observations of 3 program stars 
with and without the filter, and its extinction was mea- 
sured at 4.753±0.039 mag at J, 4.424±0.033 mag at H, and 
4.197±0.024 mag at Ks- 
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HD 129333 with PALAO/PHARO, and for HD 49197, 
V522 Per, and RX J0329. 1+0118 with NIRC2 (Matthews 
et al., in prep.) and the Keck II AO system (diffraction 
hmit 0.05" at Ks; Wizinowich et al. 2000). Conditions 
were not photometric during foUow-up, and only the best 
images (Strehl ratio S > 40%) were selected for astrom- 
etry. HD 129333 was followed up in the narrow-band 
Brackett 7 filter, which allowed higher throughput than 
the ND 1% filter in the shortest (1.4 sec) PHARO expo- 
sures. Given the unstable atmospheric conditions during 
the second-epoch imaging of HD 129333, this allowed 
us to take high signal-to-noise {S/N) exposures on time- 
scales that would most finely sample the variations in the 
seeing, and to select only the ones with the best imaging 
quality. 

Keck follow-up is done at JHKs, or at KgL' if the can- 
didate companion is expected to be bright enough to be 
seen at L' (giv en 0.7 < Ks — L' < 2.5 for L and T dwarfs; 
iLeggett et al.l 12002). Sequences of short (up to 20 sec 
from multiple co-adds) dithered non-coronagraphic, and 
long (1 min) target-sky-target exposures are taken with 
a 1"- or 2"-diameter coronagraph in the same manner as 
with PHARO, though without detector rotations. The 
candidate companions are exposed until a S/N ratio 
comparable to that in the first-epoch PHARO obser- 
vation is achieved (for similar positional accuracy), up 
to 6 min per filter in J, H, and Ks- The total in- 
tegration time at L' is up to 10.5 min, which allows 
the detection of L' < 15.0 objects. We mainly use 
the 40 mas/pix (wide) NIRC2 camera (41" FOV), which 
severely under-samples the Keck AO point-spread func- 
tion (PSF), but is known to suffer from less distortion 
than the 20 mas/pix (medium) camera over the same 
field ijThompson. Egami. &: Sawickill200H) . Although we 
also have the option of using the 10 mas/pix (narrow) 
camera (10" FOV) in NIRC2, it does not allow follow-up 
of distant candidate companions, and we avoid using it 
for consistency with the other NIRC2 observations. 

2.1.3. Imaging Data Reduction 

All imaging data are reduced in a standard fashion for 
near IR observations. Flat fields are constructed either 
from images of the twilight sky (for the Palomar data), or 
from images of the lamp-illuminated dome interior (for 
the Keck data). A bad pixel mask is created from the 
individual flats, based on the response of each pixel to 
varying fiux levels. Pixels whose gain deviates by more 
than 5 sigma from the mode gain of the array are flagged 
as bad. Sky frames for the dithered, non-coronagraphic 
exposures are obtained by median-combining four of the 
five exposures in the dither pattern (excluding the central 
pointing), and rejecting the highest pixel value in the 
stack. The coronagraphic-mode sky frames are median- 
combined using an average sigma clipping algorithm to 
remove pixels deviant by more than 5 sigma. The sky- 
subtracted images of each target are divided by the flat 
fleld, then registered, and median-combined to create a 
final high signal-to-noise {S/N) image (Figures d and EJ. 
However, photometric and astrometric measurements are 
performed on the individual reduced images. 

No PSF stars are observed at either Palomar or Keck. 
With PHARO at Palomar, median-combined images 
from all 4 detector orientations can be used to repro- 
duce an approximate PSF. This approach was chosen to 



emulate the observation of separate PSF stars of iden- 
tical brightness and color, while optimizing the time 
spent on science targets. However, we have found that 
a simple 180° rotation and subtraction technique works 
equally well, and we use that on both the Palomar 
and Keck data. While neither approach eliminates tele- 
scopic speckle noise (as could be the case if actual PSF 
stars were observed), both significantly reduce point- 
symmetric structure in the PSF. 

2.2. Astrometric Calibration 

The exact pixel scale of the 25 mas PHARO cam- 
era was determined using known binary stars from the 
Sixth Orbit Catalog faartkop f. Mason, fc Worlevll2001l: 
Hartkonf fc Ma,sonll2003al): W DS 09006+4147 (grade 1; 
Hartkopf. Mason, fc McAlister..l99a) . WDS 16147+3352 
fsjrade 4: |gca,rdial 11 9791) . WDS 18055+0230 (gra de 1: 
I Pourbaixi liMfrand WDS 20467+1607 (grade 4; Ih^ 
Il994j) . These "calibration binaries" are observed 
throughout our campaign at Palomar at all four de- 
tector orientations. The combination of grade 1 (ac- 
curately determined, short-period) orbits and grade 4 
(less accurately known, longer-period) orbits was selected 
from the list of astrometric calibrators recommended by 
iHartkopf fc Ma son (2003b). Despite the lower quality of 
the solution for binaries with grade 4 orbits, their peri- 
ods are generally much longer (889 and 3249 years for 
WDS 16147+3352 and WDS 20467+1607, vs. 21.78 and 
88.38 years for WDS 09006+4147 and WDS 18055+0230, 
respectively), so their motions are predicted with suffi- 
cient accuracy for many years into the future. The mean 
pixel scale of PHARO was measured to be 25.22 mas/pix 
with a 1 sigma scatter of 0.11 mas/pix among mea- 
surements of the individual binaries at different detec- 
tor orientations. This measurement is consistent with, 
though less precise th an our previo us determination 
(25 . 1 68 ± 0. 034 mas /pix: iMetchev. Hill enbrand. & Whiti 
I2OOI . which was obtained from only one calibration bi- 
nary at a single Cassegrain ring orientation. The larger 
scatter of our more recent measurement is indicative of 
the systematics involved in choosing different calibration 
binaries, and in observing at more than one detector ori- 
entation. 

The plate scale of the wide NIRC2 camera 
was calib rated using the binary WDS 15360+3948 
(grade 1; iSoderhieTinl I1999H . The obtained value of 
39.82+0.25 mas/pix is consistent with the pre-ship mea- 
surement of 39.7+0.5 max/pix^. Because only one binary 
was used for calibrating NIRC2, a 0.63% error term cor- 
responding to the uncertainty in the semi-major axis of 
the binary has been added in quadrature to the uncer- 
tainty of our measurement. 

2.3. Spectroscopy 

2.3.1. Observations 

Spectroscopy of faint targets in the halos of bright ob- 
jects is more challenging than spectroscopy of isolated 
targets. Care needs to be taken to achieve optimum sup- 
pression of the scattered light of the bright stars either 
by using a coronagraph (in a manner similar to corona- 
graphic imaging) , or by aligning the slit so as to minimize 

^ See |htt p : / / www2 . keck. hawaii. edu / inst /mrc2/genspecs . html | 
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light admitted from the halo. Since neither PHARO nor 
NIRC2 allow coronagraphic spectroscopy, when taking 
spectra of faint companions we orient the slit as close as 
possible to 90° from the primary-companion axis. Spec- 
tra of brighter companions, for which the signal from the 
halo is negligible compared to that from the target, are 
obtained by aligning the binary along the slit. Given that 
we often use the F5~K5 IV-V primaries in our sample as 
telluric standards, such alignment improves our observ- 
ing efficiency. Sky spectra are obtained simultaneously 
with the target spectra, by dithering the targets along 
the slit. 

Promising candidate companions are observed at 
medium resolution {R = 1000 - 3000) at K and (AO 
correction permitting) at J. In PHARO we use the cor- 
responding grism and filter combination {K or J) and the 
40 mas/pix camera, while in NIRC2 we use the K oi J fil- 
ter with the "medres" grism, and the wide (40 mas/pix) 
camera. For ii'-band spectroscopy in first order with 
PHARO and in fifth order with NIRC2 we achieve com- 
plete coverage of the 2.2fim atmospheric window. For 
J-band spectroscopy in first order with PHARO the cov- 
erage is limited by the size of the detector to a 0.16^m 
bandwidth. ,/-band spectroscopy with NIRC2 can cover 
the entire I.2/im window, but the spectrum is split be- 
tween the fifth and the sixth dispersion orders. To allow 
simultaneous data acquisition from both orders, we fit 
both onto the wide camera by clipping the fifth (longer- 
wavelength) order shortward of 1.22/im and the sixth 
(shorter-wavelength) order longward of 1.28/im. The 
combined J-band spectrum has complete coverage be- 
tween I.16-1.35/im. 

We took J and if-band spectra with PHARO 
of HD 129333A/B, J and K spectra with NIRC2 
of V522 PerA/B, and X-band NIRC2 spectra of 
HD 49I97A/B and RX J0329.I+0118A/B. For the 
observations of HD 129333 with PHARO, and of 
RX J0329.1-hOII8 with NIRC2, the binary was ahgned 
along the slit, whereas spectra of the other two binaries 
(with fainter companions) were obtained by placing each 
component in the slit individually. A 0.26" (6.5 pix) slit 
was used in PHARO, and a 0.08" (2 pix) slit was used 
in NIRC2, resulting in spectral resolutions of 1200 at J 
and 1000 at K with PHARO, and 2400 at fifth-order J 
(1.22-1.35^m), 2900 at sixth-order J (1.15-1.28/im), and 
2700 at K with NIRC2. 

Flat fields with the dispersive grisms in place were not 
obtained for any of our spectroscopic observations. In- 
stead, the raw spectra were divided by an imaging flat 
field, constructed in the same manner as for the imaging 
observations. The flat-fielded spectra were corrected for 
bad pixels, using the same bad pixel mask as in the imag- 
ing case. Strong positive deviations due to cosmic ray 
hits were then elimin ated using the L.A. Cosmic Lapla- 
cian fllter algorithm l|van Dokkumll200l . Fringing was 
at a noticeable level (^15%) only in spectra obtained 
with PHARO (of HD 129333), and was reduced to below 
5% by dividing the target spectrum by that of the telluric 
standard, taken at a nearby position on the detector. 

2.3.2. Spectroscopic Data Reduction 

Extraction of the spectra is performed using 
iraf/apextract tasks. To reduce contamination from 
the halo of the primary in the companion spectra, the 



local background (arising mostly from the halo) is flt by 
a flrst-order polynomial along the direction perpendicu- 
lar to the dispersion axis and subtracted during extrac- 
tion. In addition, the aperture width for the companion 
spectrum is conservatively set to the full width at half 
maximum of the proflle so as to include only pixels with 
maximum signal-to-noise. Pixels within the aperture 
are summed along detector columns (which are nearly 
perpendicular to the dispersion axis), and the resulting 
"compressed" spectrum is traced along the dispersion 
axis by fitting a high order (8 to 15) Legendre polyno- 
mial. The tracing step is approximately equal to the slit 
width, where at each step the data point is obtained as 
the sum of 3-10 adjacent pixels (1.5-3 slit widths; the 
spectra of fainter objects being more heavily averaged) 
in the compressed spectrum. This procedure is found 
to produce consistent results for extractions of multi- 
ple dithered spectra of the same companion, indicating 
that scattered light contamination from the primary has 
been reduced to a small level. Nevertheless, in some 
cases of very faint close-in companions, the continuum 
shape is still found to vary noticeably among the individ- 
ual extractions. We therefore avoid classifying the spec- 
tra of the companions based on their continuum shapes, 
but rely on the relative strengths of narrow absorption 
features (discernible given our resolution) instead (Sec- 
tion ESI- 

A dispersion solution for each spectrum is obtained 
from night sky lines in non-sky-subtracted images, using 
the task identify and OH emission-line lists available 
in IRAF. For the spectra of primaries observed sepa- 
rately (e.g., HD 49197 and V522 Per), whose shallower 
(10-120 sec) exposures do not contain telluric emission 
lines at a high enough S/N to allow the fitting of a dis- 
persion solution, such is derived after registration with 
deeper (10-15 min long) companion spectra taken im- 
mediately after those of the primaries. The tasks FXCOR 
and SPECSHIFT are used to cross-correlate and align the 
individual wavelength-calibrated spectra for a given ob- 
ject. The individual primary and companion spectra are 
then median-combined using SCOMBINE. 

Since the primaries (when earlier than KO) double as 
telluric standards, their spectra are first corrected for 
photospheric absorption from atomic hydrogen (H P/3 at 
1.282/im and H Br7 at 2.166/J,m) by hand-interpolating 
over the absorption with the task SPLOT. Weaker ab- 
sorption features due to Na I, Ca I, Si I, Al I, and Fe I, 
although present in the spectra of our telluric standards, 
are left uncorrected because of blending (at our resolu- 
tion) with various OH absorption lines in the Earth's 
atmosphere. After thus correcting the telluric spectra, 
the spectra of the companions are divided by those of 
the telluric standards. Finally, the spectra of the com- 
panions are multiplied by black bodies with temperatures 
corresponding to the spectral types of the telluric stan- 
dards (ba sed on effe ctive temperature vs. spectral type 
data from Coxl200(il) , and boxcar-smoothed by the width 
of the slit. 

The above reduction procedure applies in exact form 
for NIRC2 spectra, which do not suffer from fringing. For 
PHARO data, which are noticeably affected by fringing, 
we divide the individual target spectra by those of the tel- 
luric standard before median-combining and wavelength 
calibration, since fringing depends on detector position, 
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not on the dispersive element. 

Reduced spectra for the objects reported here are pre- 
sented in Figure |2| and 01 The mismatch in the contin- 
uum slopes between the fifth- and sixth-order the J-band 
spectra of the companion to V522 Per may be due to our 
use of imaging, instead of spectroscopic fiats. 

3. ANALYSIS 
3.1. Photometry 

Broad-band near-IR photometry of the companions is 
presented in Table|3| The measured quantity in each case 
is the relative flux (AJ, AH, AKs, AL') of the com- 
panion with respect to that of the primary. When the 
companions were visible without the coronagraph (in all 
cases except for the close-in companion to HD 49197), the 
fluxes of both components were measured directly from 
the short-exposure non-coronagraphic images. Flux ra- 
tios for the HD 129333, V522 Per, RX J0329.1-I-0118, 
and HD 49197A/"C" systems were obtained from the 
Palomar images in apertures of radii of 0.2" {2X/D at 
Ks). For the photometry of HD 129333B we subtracted 
the halo of the nearby (0.74") primary (as detailed in 
Section I2.1.3|l to minimize its contribution to the flux 
of the secondary. In the case of the HD 49197A/B sys- 
tem, the 0.95" companion is not seen in PHARO images 
taken without the coronagraph, and only coronagraphic 
exposures were obtained with NIRC2. The magnitude of 
the companion was obtained with respect to the residual 
flux of the primary seen through the 1" NIRC2 coron- 
agraphic spot. The flux measurements were performed 
on the PSF-subtracted NIRC2 images in a 0.16"-radius 
{3.2X/D at Ks) aperture. The transmission of the spot 
was established from photometry of another program 
star in images taken with and without the 2" corona- 
graph: its extinction was measured at 9.27±0.07 mag at 
J, 7.84±0.03 mag at H, and 7.19±0.03 mag at Ks. The 
companion was seen in all 6 one-minute coronagraphic 
exposures at Ks, but because of varying photometric 
conditions and proximity of diffraction spikes, it was de- 
tected in only 3 out of 6 PSF-subtracted exposures at H, 
and 1 out of 6 at J. 

Sky values for each of the observed objects were deter- 
mined as the centroid of the flux distribution in a 0.1"- 
wide annulus with an inner radius larger by X/D than the 
distance at which the radial profile of the object fell be- 
low the level of the local background. For the primaries 
observed with PHARO this inner radius was 2.0" (20A/I? 
at Ks), while for HD 49 197 A, the flux of which was mea- 
sured through the NIRC2 coronagraph, the inner radius 
was 0.2" {5X/D). For the companions, the inner radius of 
the sky annulus varied from 2.25-4 times the Ks diffrac- 
tion limit, depending on whether the local background 
was strongly influenced by the halo of the primary (as 
near HD 49197 and HD 129333), or not. Experiments 
with varying sky annulus sizes in the two more detached 
systems (V522 Per and RX J0329.H-0118) showed that 
the relative photometry between two objects in the same 
image is preserved to within 0.08 mag for annuli ranging 
between 2.25-20A/D in inner radius. 

The J — H and H — Ks colors, and Ks magnitudes, are 
derived from the measured relative photomet ry by adopt- 
ing th e Two-Micron All Sky Survey f2MASS: iCutri et all 
|2003^ magnitudes for the primaries. The Ks — L' color 



of the companion to V522 Per was calc ulated assuming 
Ks - L' = 0.04 for the F5 V primary l|Bessell fc Bretl) 
IT988I) . A near-IR color-color diagram of the detected 
companions is presented in Figure |S| Extinction cor- 
rections of Aj = 0.087, Ah = 0.055, Ak = 0.033, and 
Al = 0.016 have been applied to the colors of V522 PerB, 
based on E{B — V) = 0.10 toward the a Per cluster 
as measured by iPinsonneault et aTl l)1998|) . We have 
adopted A y = 3.1E(B — V), and the in terstellar extinc- 
tion law of lCardelh. Clavton. fc Math is (1989'). 

Based on the colors and the AKs magnitudes, we can 
infer that HD 49197 "C" is a likely F-G background star, 
whereas the remaining companions are consistent with 
being late- type stars associated with their primaries. 

3.2. Astrometry 

Three of the four systems were observed at two dif- 
ferent astrometric epochs. RX J0329. 1-1-0118 was ob- 
served only onc e, though prior-epo ch astrometric data 
for it exist from lSterzik et all l|1997D . Table El details for 
each binary the observed offset and position angle of the 
companion from the primary during the first and second 
epochs of observation, as well as the estimated offset and 
position angle during the first epoch, had the system not 
been a common proper motion pair. The estimates are 
extrapolated backwards in time from the second-epoch 
astrometry, which for all except HD 129333 was ob- 
tained with NIRC2 on Keck and was m ore accurate, as- 
sumin g proper motions from Hipparcos (|Perrvman et all 
E997) for HD 49197 and HD 129333, and from Tycho 2 
lH0g et al. 2000) for V522 Per and RX J0329. 1+0118. 
Parallactic motions were also taken into account, as they 
are significant for stars < 200 pc at our astrometric pre- 
cision (several milli-arcsec). The assumed proper mo- 
tions and parallaxes for all stars are listed in Table |2| 
To ensure the correct propagation of astrometric er- 
rors, the epoch transformations were performed following 
the co-var i ant tr eatment of the problem, as detailed in 
lLindegrenl l )1997j) . 

Proper motion diagrams for each object are presented 
in Figures EHHl The first- and second-epoch measure- 
ments are shown as solid points with error bars, and 
the inferred first-epoch position (assuming non-common 
proper motion) is shown only with errorbars. The dashed 
line reflects the combined proper and parallactic motion 
of the primary over the period between the two epochs. 
Below we discuss the evidence for common proper motion 
in each system. 

3.2.1. HD 49197 

The existence of the close (0.95") companion to 
HD 49197 was unappreciated prior to the second-epoch 
imaging: the star was followed up because of the more 
distant (6.8") candidate companion (HD 49197"C"). 
Upon its discovery in the Keck image (Figure H^), the 
close companion was recovered in the earlier Palomar 
image, where a dark circular ring around the core at the 
radial distance of the first Airy null («0.1") distinguishes 
it from the telescopic speckles (Figure Qd). 

As is evident from FigureEl HD 49197B (left panel) is 
much more consistent with being a proper motion com- 
panion of the primary, than with being an unrelated 
background star projected along the same line of sight, 
whereas the reverse holds for HD 49197"C" (right panel). 
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Therefore, we claim that the close-in HD 49197B is a 
bona-fide companion, whereas the more widely-separated 
HD 49197"C" is not. 

The astrometry for HD 49197B from the two ob- 
servational epochs is not in perfect agreement, per- 
haps because of its orbital motion around component 
A. At a projected separation of 43 AU, the orbital pe- 
riod will be >240 years (assuming a circular face-on 
orbit and a mass of 1.16M(7i f or the F5 V primary; 
lAllende Prieto &: Lambe rr'1999'1. resulting in a change 
in position angle of <2.6° between the two epochs. This 
may explain the observed discrepancy of 0.64° ± 0.47° in 
position angle between the two observations, while the 
change in separation (2.4±5.8 mas) is consistent with 
being zero. Future observations spanning a sufficiently 
long time-line may be used to determine the dynamical 
mass of the HD 49197A/B system. 

3.2.2. HD 129333 

The offset positions of the companion to HD 129333 
in the two observations taken 16 months apart are fully 
consistent with each other, and inconsistent with the ob- 
ject being a background star (Figure |7|. HD 129333B 
thus shares the proper motion of the primary, and is a 
bona-fide companion. 

3.2.3. V522 Per and RX J0329. 1+0118 

The astrometry for these two systems (Figure IH)) is in- 
conclusive because of their smaller proper motions, insuf- 
ficient time-span between the observations, and/or less 
accurate astrometry. The likelihood of physical associa- 
tion of the companions with the primaries is investigated 
from low-resolution spectroscopy in Section l4. li b elow . 

3.3. Spectroscopy 

Infrared spectral classification of M-L dwarfs is done 
most successfully in the J and H bands, where a 
suite of indices based on the relative strengths of 
H2O, FeH, K I, and Na I absorption have been 
developed to characterize their effe ctive temperatures 
jSlcsnick. Hillenbrand, fc Carpcnter"2004': 'Gcball e et all 
12002; .Gorlova et al...2003t iMcLean et al...200a) and sur- 
face gravities l|Gorlova et al.l I2003D . However, spec- 
troscopy of cool companions in the bright halos of their 
primaries is often more difficult at J and H than at 
if, because of larger flux contrast and poorer quality 
of the AO correction at shorter wavelengths. Here we 
present J-band spectra of the two brighter compan- 
ions: HD 129333B, and V522 PerB. Cool dwarfs can 
be classified at K band based on the strength of their 
H2O absorption shortward of 2.05/im, the depth of the 
CO 2.29/im bandhead, and the equivalent width (EW) 
of the Na I 2.21/xm double t (Klcin mann fc Hal|ll986t 
lAh et alJlT995t IMcLean et al.ii2003tiSlesnick et al.ll2004H . 
Even cooler, late-L and T dwarfs, are best character- 
ized by the stre ngth of t heir CH4 absorption at >2.20/im 
(|Burgasser et~al.„.2002: .McLean et aLll2003D . We present 
if-band spectra of all four companions discussed in this 
paper. 

Spectral types of the detected companions have been 
determined following infrared absorption line classifi- 
catio n systems at i^T-ba nd (for Na I, Ca I, CO) 
from IMcLean et al.l l)2003l spectral types M6V-T8) and 



lAli et alJ |IT99RI F3V-M6V), and at J -band (for K I 
and Ti I) fr o m iGorlova et alJ l|2003[ M4V-L8) and 
I Wallace et all l)2000l 07-M6'). Luminosity classes are 
based on the relativ e strengths of Na I, Ca I, and CO ab- 
sorption at K band ( Kle inmann fc Hall 1986), and on the 
stren gths of Mn I, an d K I absorption at J ( Wallace et alJ 
I2ti00t IGorlova et al.ll2f)n^ . 

We have avoided the use of temperature-sensitive H2O 
indices that span large fractions of the s pectrum (>4% 
total band width, e.g., the H2OD index of McLea n et alJ 
((2003i) . or the H2O-2 index of Slc snjck et al, (2004')) be- 
cause of their strong dependence on the overall con- 
tinuum shape of the spectrum, and because of the 
uncertainties in the spectral shapes of faint compan- 
ions extracted from the halos of bright objects (Sec- 
ti on 12.311. We have used the J-band 1.34 /im H2O index 
of IGorlova et al.l ()2003j) however, which measures only 
the onset of H2O absorption at that wavelength, and is 
narrow (1.4% band width). 

Below we analyze the K- and J-band spectra of the 
companions to HD 49197, HD 129333, V522 Per, and 
RX J0329. 1-1-0118. The inferred effective temperatures 
and spectral types are provided in Tables and 

3.3.1. K-hand Spectroscopy of HD 49197B 

Our spectrum of HD 49197B (Figure 13 shows strong 
CO and H2O bands characteristic of ultra-cool dwarfs, 
but lacks a CH4 absorption feature, indicating that it is 
earlier than L8 (Gcballe ct al. 2002). On the other hand, 
Na I absorption is also absent from the spectrum, p oint - 
ing to a spectral type of L3 or later (IMcLean et al.l2003..'l 
indep endent of gravity (cf. Fig^ire 8 in IGorlova et all 
I2003|) . Following the analysis of i?-2000 if -band spec- 
tra of M6-T8 dwarfs in McLean ct al. (2003), we form 
a CO absorption index from the ratio of the median 
flux in an absorbed (2.298-2.302 /im) vs. an unabsorbed 
(2.283-2.287 /im) region of the spectrum. We find 
CO=0.80±0.03 (where the error has been estimated as 
the quadrature sum of the relative uncertainties of the 
median in the two spectral regions), which indicates a 
spectral type of L7 or later. However, the CO index is 
not very sensitive to L-dwarf temperatures, and varies by 
up to 0.15 (60% of its total range of variation between 
spectral type s M5 and T2) within the same spectral type 
(Figure 13 in IMcLean et al.ll2003t) . 

Alternatively, we can also use the absolute A'5-band 
magnitude of HD 49197B to estimate its spectral type, 
follow ing the empirical relation of iKirknatrick et afl 
l20tir)j) ■ based on a sample of 24 M and L dwarfs with 
measured parallaxes: 

Mks = 10.450-l-0.127(subclass)-h0.023(subclass)^ (1) 

where subclass = -1 for M9 V, for LO V, 1 for LI V, 
etc. The scatter about the fit is approximately 1 sub- 
class. From the inferred apparent Ks magnitude of 
HD 49197B, and from the parallax of the primary, we 
obtain = 11.04 ± 0.24 for the secondary, which cor- 
responds to a spectral type of L3±1.5. 

We assign a final spectral type of L4 with an uncer- 
tainty of 1 subclass. This is based on the intersection 
of the results from our spectroscopic analysis, suggest- 
ing L3-L7, and from the ifg-band absolute magnitude, 
pointing to L1.5-L4.5. A spectral type of L4±l for 
HD 49197B is also consistent with a by-eye comparison of 
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its K-hand spectrum with the grid of L-dwarf standards 
from lLeggett et alJ 1)20011) . 

3.3.2. K-band Spectroscopy of HD 129333B, 
RX J0329.1+0118B, and V522 PerB 

We classify the if-band spec tra of th e se com pan- 
ions following the anal yses of lAli et all 1)1995)) and 
iKleinmann fc Halll l)1986)l . whose data are taken at simi- 
lar resolutions to ours (i?=1380-3900 and 2500-3100, re- 
spectively) , and span the F8-M7 spectral type range. We 
employ the spectral classification sequence of lAli et al.l 
whose empirically-derived in dices are based on a larger 
sample of dwarf stars than in IKleinmann fc Hall 

Our reduced spectra were first shifted to km s~^ 
heliocentric velocity, where the shift was determined by 
fitting Gaussian profiles to the Na I doublet, and com- 
paring the fit ted line center s to their values in the so- 
lar spectrum l|MohleillT955)) . For consistency with the 
lAli et al] l)1995[l spectral classification, we have chosen 
the same bands for integrating the Na I (2.21/im) Ca I 
(2.26^m), and 1^00(2-0) (2.29/im) absorption. The con- 
tinuum in the target spectra was fit to three regions 
devoid of absorption lines: 2.0907-2. 0951^m, 2.2140- 
2.2200/.tm, and 2. 2873-2. 2925/xm. These have been se- 
lecte d as a c ombination of the continuum regions used 
bv .Ali et al.l l(T995^) and J^leinmann & Hah (1986), so as 
to constrain the fit on both sides of the Na I doublet (as 
in IKleinmann fc Hal]|ll98'^ where the continuum-fitting 
regions are wide ly-separated fro m the Na I lines) , as well 
as near it (as in'A li et al.lll995[ where the continuum is 
constrained only on the long-wavelength side of the Na I 
lines) . 

The absorption strength in each band was obtained 
as an equivalent width, by integrating the profile of the 
spectrum in the band with respect to a global contin- 
uum level defined by the three continuum bands. The 
only exception is the EW of the CO (2-0) first over- 
tone bandhead, for which we have adopted the mean 
continuum level of the thir d continuum b and (as in 
IKleinmann fc Hal]l ll98ff and lAli et al.lll995D . The one- 
sigma errors in the EWs were calculated by propagating 
the r.m.s. noise of the spectrum in the nearest continuum 
band, assuming independent pixel variances. The EWs 
of Na I and Ca I were corrected for corresponding ab- 
sorption in the spectra of the telluric standards, the EW 
of which (0 84±0.57 A for Na I and 1.14±0.36 A for Ca I 
in G stars; lAli et aLlll995)l was added to that measured 
for the companions. The final EWs are listed in TableEl 

To infer effective temperatures for the companions to 
HD 129333, RX J0329. 1-^0118, and V522 Per, we em- 
ploy empirical relations between the strength of K-hand 
Na I, Ca I, and CO absorptio n and ef f ective tempera- 
ture (Teff), as determined bv lAli et al.l l)1995l their Ta- 
ble 4). Na and CO produce the most characteristic 
K-haiid features of cool stars. However, their absorp- 
tion strengths are inaccurate tracers of temperature for 
stars cool er than 4000 K, and are , in addition, gravity - 
sensitive l)Kleinmann fc Hal]llT98^ iGorlova et al.ll2003l) . 
On the other hand, Ca I transitions in the if -band re- 
quire the population of higher energy states, and hence 
are more temperature-sensitive. However, their absorp- 
tion strength is degenerate with the stellar effective tem- 
perature: Ca I absorption at 2.26/im peaks at ~3500 K, 
and decreases for higher, and lower effective tempera- 



tures l)Kleinmann fc Hal]llT98^IAn et al.llT99l . This be- 
h ayior is fit via t\y o separate Teg vs. EW(Ca I) relations 
in lAli et al.l l)1995)) . By combining the information from 
Ca I with that from Na I and CO (2-0), we can break 
this degeneracy, and use the more temperature-sensitive 
Ca I index to constrain the effective temperature for each 
star to w ithin ~300 K (t he quoted uncertainty of the Ca I 
index in lAh et alJll995)) . For HD 129333B, with a Ca I 
absorption strength near the breaking point between the 
"hot" and "cool" relations, we take the average of the 
two estimates. For V522 PerB and RX J0329.1-^0118B, 
whose Na I and CO absorption is indicative of temper- 
atures Toff < 3400 K, we use only the "cool" relation. 
The effective temperatures inferre d from C a I ab sorp- 
tion are listed in Table [S] As in lAli et all l)1995)) , we 
adopt a spectral type vs. effective temperature classi- 
fication for M dwarfs from Bessell (1991), and obtain 
spectral types of Ml, M3, and M4 for HD 129333B, 
RX J0329.1-h0118B and V522 PerB, respectively, with 
an uncertainty of 2 spectral subtypes. A visual inspec- 
tion and comparison of the strengths of the various ab- 
sorption featu res with if -band spectral seq uences from 
iLeggett et all ()2001)) and iWallace fc Hinklei l)1997t) con- 
firms these results. Given the comparable values of their 
Na I and 1^00(2-0) EWs, aU stars have likely dwarf 
gravities (cf. Figure 7 in Kleinmann & Hall 1986). 

3.3.3. J -band Spectroscopy of HD 129333B and 
V522 PerB 

We further constrain the spectral types of HD 129333B 
and of V522 PerB from their J-band spectra. These 
show the characteristic absorption features of M stars: 
K I lines at 1.169/.tm, 1.177/xm, 1.243/im, and 1.252/im, 
Fe I at 1.188/im and 1.197^m, Na I at 1.268^m, Al I 
at 1.312/im and 1.315/im, and H2O absorption at A > 
1.34^m (Figure0J. The spectrum of V522 PerB exhibits 
also Ti I and Mn I absorption over 1.282-1.298 /im. 
Spectral classification at J band was done based on 
the depth of the H2O and K I absorption , following 
the index definitions of IGorlova et al.l l)2003f) . After re- 
sampling oui_dx]ia^_toniatch the i?«350 spectral reso- 
lution of IGorlova et all we form the 1.34/im water in- 
dex as the ratio of the mean fluxes (in a 0.004 fxm wide 
region) at 1.336/xm and at 1.322/im, and we measure 
the K I EW by integrating the absorption over the re- 
gion 1. 2375-1. 2575/im. Because K I absorption in the 
solar-like photospheres of the telluric standards is small 
{EWiK I A1.14Mm) ^ 0.1x gVF(N a I A2.21^m) = 0.08 A 
for the Sun; MohleJ IT955t lAH et al. 1995), the K I EW 
measurements of the companions were not corrected for 
it. 

For HD 129333B, we cannot estimate the strength 
of the 1.34/im water absorption because of insufficient 
spectral coverage. The EW of K I indicates a spectral 
type of M2-M4. Averaging this with our K-hand es- 
timate of Ml±2, we assign a spectral type of M2±l for 
HD 129333B. Given the youth of HD 129333 (<120 Myr; 
see Section 14.211 , the companion may have lower-than- 
dwarf gravity. The effect of this on alkali absorption 
lines in the n ear IR is degenerate with temperature 
l)Gorlova et all 12003: McGovern et al. 2004), and could 
be compensated by a later spectral type. However, spec- 
tral types later than M3 are inconsistent with the depth 
of the strongly temperature-sensitive Ca I absorption in 



8 



Metchev & Hillenbrand 



this star fSection l3.3.2|) . In addition, HD 129333B lacks 
noticeable Mn I absorption, which is weak in M dwarfs 
but grows deeper with decreasing surface gravity in M 
stars (Fig. 9 in Wallace et al. 2000). Hence, we conclude 
that HD 129333B is a M2±l dwarf. 

For V522 PerB, both the K I EW and the water absorp- 
tion index point to a spectral type of M3-M5, in agree- 
ment with our ii'-band estimate (M4±2). The stronger 
Ti I absorption than in HD 129333B is also consistent 
with a cooler photosphere. We thus assign a spectral 
type of M4±l to V522 PerB. The strength of the Mn I 
transition indicates a potentially subgiant surface grav- 
ity, though as noted in Section 13.3.21 the effect is not 
seen at K band. 

4. DISCUSSION 

4.1. Likelihood of Physical Association 

Our astrometric follow-up of HD 49197B and 
HD 129333B confirmed common proper motion between 
these two companions and their respective primaries. 
However, the smaller proper motions of V522 Per and 
RX J0329. 1-1-0118 prevented us from concluding the 
same for their respective companions, given the time- 
span of our observations. The probability of physical 
association in these systems can be inferred from the 
spectroscopically determined spectral types and absolute 
magnitudes of the companions. If the absolute mag- 
nitude inferred from the spectral type of a companion 
agrees with its measured apparent magnitude at the he- 
liocentric distance of the respective primary, then the 
companion is likely to be a bona-fide one (modulo the 
space density of stars of similar spectral type as the com- 
panion) . 

Figure El presents a comparison of the spectroscopic vs. 
photometric absolute magnitudes. The correspondence 
is good for the four companions followed up via near IR 
spectroscopy, indicating that they are at similar helio- 
centric distances as their primaries, and are thus likely 
to be physically bound to them. The location of the re- 
maining candidate companion (HD 49197"C") along the 
ordinate is inferred from its near-IR colors. As men- 
tioned in Section EH HD 49197"C" is a likely F-G star 
(2.0 < Mk < 4.0; Cox 2000), i.e., it is intrinsically too 
bright to be associated with HD 49197 (F5 V) given its 
faint apparent magnitude. 

A robust statistical analysis of the likelihood of chance 
alignment in the four systems discussed here is not yet 
possible at this stage of the survey. They are only a frac- 
tion of the ones discovered to have candidate compan- 
ions. With follow-up observations still in progress, the 
exact number of bound systems is unknown. We defer 
a discussion of the companion chance alignment proba- 
bility based on the full ensemble statistics until a later 
paper. Here we consider these probabilities only on a per 
system basis. To give an approximate idea of the lim- 
ited statistics from which these preliminary results are 
extracted, we point out that to date we have analyzed 
multi-epoch astrometric data for approximately 15 stars 
(mostly from the deep survey) with faint (AKs > 3, i.e., 
expected to be of spectral type M or later) companions 
within 4". 

We base our calculation of the probability of false as- 
sociation in each system on the empirically determined 



spatial density of cool objects (M-T spectral types) in 
the solar neighborhood. There are 112 such known ob- 
jects in the northern {6 > —20°) 8 pc sample (jReid et alJ 
The northern 8 pc sample covers 65% of the sky, 
and is estimated to be ^15% incomplete. The total num- 
ber of cool objects and white dwarfs within 8 pc of the 
Sun is therefore expected to be 198, with a volume den- 
sity of 0.10 pc~'^. This estimate is based on a small 
fraction of the thin disk population (scale height 325 pc; 
iBahcall fc ^Soneira 1980) of the Galaxy and hence should 
not vary substantially as a function of galactic latitude. 

We then calculate the number of cool dwarfs expected 
to be seen in projection toward each system within a con- 
ical volume of radius 4" centered on the star, with the ob- 
server at the apex of the cone. We truncate the radial ex- 
tent of the conical volume by requiring that the apparent 
K magnitude of a projected companion falls within the 
limits allowed by the spectral type (and hence, absolute 
magnitude) of the detected one. Absolute K magnitudes 
for the M2-4 dwarfs discussed here have been adopted 
from Bessell (1991). Although Bessell's M dwarf classifi- 
cation system pre-dates the discovery of ultra-cool dwarfs 
(later than M5), it remains valid for early M dwarfs. For 
L4± l spectral types w e adopt absolute K magnitudes 
from iDahn et all l)2002() . 

The expected number /i of unrelated cool dwarfs within 
the relevant volume around each star is listed in the last 
column of Table El Given that for all stars ^ ^ 1, we 
can assume that the event of seeing an unrelated field 
object in the vicinity of one of our program stars is gov- 
erned by Poisson statistics. Hence, the probability of 
finding one or more such dwarfs near any given star (i.e., 
the probability of chance alignment) is 1 — ~ /j,. As 
seen from Table after having followed up the compan- 
ions spectroscopically, we can claim with >99.99% cer- 
tainty in each case that the companion is physically as- 
sociated with its respective primary. As discussed above, 
such probabilities need to be regarded in the context of 
the ensemble statistics. Within our sample of 4 spectro- 
scopically confirmed companions, the probability that at 
least one is a false positive is 3 x 10"'^. This exemplifies 
the power of spectroscopic follow-up as an alternative to 
multi-epoch astrometry in constraining the likelihood of 
physical association in a system. 

4.2. Stellar Ages and Companion Masses 

We estimate the ages of the primaries (Table j^l from 
published data on their chromospheric activity, Li I 
equivalent width, and kinematic association with young 
moving groups. Masses for each of the companions 
(Table jBJ were determined either from the low-mass 
pre-rn ain sequence evolutionary models of iBaraffe et alJ 
ill 9981). or from the brown-dwarf cooling mod els of 
iChabrier et al] pOOOl "DUSTY") and Burrows ^eTall 
{199 7!). We have not u sed the dust-free "COND" models 
of iBaraffe et al.l l|2003D . since they are more appropriate 
for temperatures <1300 K (i.e., for T dwarfs) when all 
grains are expected to have gravitationally settled below 
the photosphere. 

HD 491 97B. — From the strength of Ca H & K 
core emission in Ke ck/HIRES spectra of the primary, 
iWright et aO 1)20041) deter mine an age of 525 Myr 
for HD 49197, which we assume accurate to within 
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~50%, given the vari ation in chromosp heric activ- 
ity of solar- type stars (|Henrv et alJ 1199611 . No other 
age-related indicators exist in the literature for this 
F5 star. From our own high-resolution optical spec- 
tra, we measure EW(Li A6707.8) = 80 mA (Hillen- 
brand et al., in prep.), consiste nt with a Plei ades-like 
(120 Mvr: lStauffer. Schultz. fc Ki rkpatric kl 1998(1 or older 
age. Hence, we adopt an age of 260-790 Myr for the pri- 
mary. Assuming c o-evality, the mass of the secondary 
is_0X)60to°ollMQ (|Burrows et al.l IT997t iChabrier et all 
120001) . where the range of masses accomodates the one 
sigma error in the inferred absolute magnitude of the 
secondary, and the allowed age range for the primary. 
HD 49197B is thus a brown dwarf. 

HD 129333B. — The primary is a well-known 
young star, kinematically belonging to the Local 
Association (Pleiades moving group, 20-150 Myr; 
iSoderblom fc Clements! 1987: M ontes et al.ll2001h() . Re- 
sults from the Mount Wilson spectroscopic survey 
l[Soderbl omlll98^ and from the Keck/Lick r.v. program 
fWright et all 120041 show strong Ca II H & K emis- 
sion, indicat ing high leve ls of chromospheric activity 
and youth. iWright et alJ list an age of <10 Myr for 
this star, though the chromospheric activity-age rela- 
tion is not reliable for stars that young, in part be- 
cause of the large variance in rotation rates of stars 
younger than 50-80 Myr (e.g., Soderblo m ct.aL 1993), 
and b ecause the relation is n ot calibrated at such young 
ages. iMontes et all (|2001aD report strong Li I absorp- 
tion [EWilA I) = 198 mA), and conclude that the star 
is "significantly younger" than the Pleiades (120 Myr; 
l^tauffer ct al, 1998) ■ Assuming an age of 10-100 Myr 
for the system, we estimate the mass of the secondary at 
0.20;^;i^^MQ (from models of lBaraffe et 71111991 . 

V522 PerB. — The primary is a member of the a Per 
open cluster, confirmed by photometry, kinematics, and 
spectroscopy (lProsseilll992|) . From high-resolution spec- 
troscopy and determinati on of the lithi um d epletion 
boundary in the cluster, IStaiiffer et all l|1999t ) deter- 
mine an age of 90±10 Myr, consistent with a recent 
age estimate from u pper main-sequence turn-off fitting 
((Ventura et al.l ll998V Using the Lyon group stellar evo- 
lution models (Baraffe et al. 1998), we determine a mass 
of 0.085-0.15 Mq for the seco ndary. However, from their 
sub-stellar "DUSTY" code llChabrier et a,ljl2000ll . the 
treatment of dust opacity in which may be more appro- 
priate for this cool (^3200 K) star, we find that its mass 
is >0.10 Mq. We thus estimate 0.10-0.15 Mq for the 
mass of V522 PerB. 

RX .10329. 1+0118B. — iNmlhauser et alJ list 
RX J0329. 1+0118 (GO IV) as a fast-rotating {vsiui = 
70 km/s) X-ray source south of Taurus, with a Li I 
equivalent width of 110 mA: all indicators of rela- 
tive youth. Assuming a common origin and dis- 
tance wit h the stars in the Taurus molecular cloud 
(140 pc; iKenvon. Dobrzvcka. fc HartmannI I1994D . the 
authors claim that its bolometric luminosity is higher 
than that of a zero-age main-sequence star, and the 
star is therefore likely in the pre-main sequence (PMS) 
phase. From a proper- motion survev of PMS stars in 
Taurus- Auriga however. iFrink et alJ lll997i1 find that the 
young stars south of Taurus discussed in lNeuhauser et all 



l(1995(l are kinematically unrelated to those in Taurus, 
and that star formation in the two complexes must have 
been triggered by different events. From the Pleiades-like 
Li I equivalent width of RX J0329. 1-1-0118, we assign an 
age of sal 20 Myr for this star. Given the spectral type 
of th e secondary, its mass is O.2O1q'JqM0 l(Baraffe et alJ 
IT991 . 

4.3. HD 129333: Binary or Triple? 

The existence of a stellar companion to HD 129333 
has already been inferred in the r.v. work of 
[Duqucnnov fc Mayor (1991, DM91), who find that the 
star is a long-period single-lined spectroscopic binary 
(SBl). From their derived orbital parameters, the au- 
thors determine a secondary mass M2 > 0.37Mq, and 
suggest that the star be targeted with speckle interfer- 
ometry to attempt to resolve the companion. We should 
therefore consider whether the companion that we have 
resolved (and named HD 129333B) is identical to the 
spectroscopically inferred one. 

4.3.1. The Combined Radial Velocity and Astrometric 

Solution 

Combining r.v. and astrometric data presents a power- 
ful approach to fully constrain all the orbital elements of 
a binary system. In this Section we test the hypothesis 
that the DM91 and the imaged companions are identical 
by attempting to solve for the parameters of the relative 
orbit and checking for consistency with all available data. 

The orbital parameters that can be determined 
through r.v. monitoring of an SBl are the eccentricity 
e, the period P, the epoch Tq of periastron, the longi- 
tude u! of periastron, the systemic radial velocity Wrad,o, 
and the primary velocity semi-amplituide Ki. Ki is re- 
lated to the other orbital parameters through the mass 
function 

(Af,sin^)3 _ {l-e^/^PKl 
>^{M,+M2f^ 27rG ' 

where Mi^2 are the component masses, i is the inclina- 
tion of the orbit with r espect to the observer, and tt and 
G are constants fe.g.. iHeiiit^ 119781 p. 80). The orbital 
inclination i cannot be constrained from r.v. monitoring; 
hence the masses and the semi-major axes ai 2 of the 
binary components are degenerate with i. 

From astrometric observations we can fit for e, P, Tq, 
w, i, ai,2 (and hence, Mi_2), and for the only remain- 
ing parameter — the longitude of the ascending node 
Therefore, by performing a least-squares fit to the 
combined and appropriately-weighted r.v. and astromet- 
ric data, one can fully determine the orbit of a binary 
and att ain greater pre cision in estimating the orbital el- 
ements l(Morbevlll975ll . 

We first list the orbital parameters that have been 
already determined. Based on the r.v. measurements 
shown in Figure UUI DM91 find e = 0.665 ± 0.023, 
To = JD 2446932±20 = year 1987.37, uj = 188.0 ± 5.2°, 
Ki = 5.09 ± 0.20 km s"\ and P = 4575 days = 
12.53 years. DM91 state however, that the period is 
probably accurate only "to the nearest unit of logP" 
(i.e., lO^-^ < P < lO^-^ days, or between 8.7 and 87 
years), and calculate the uncertainties in e, Tqj and 
Ki for a fixed P. Nevertheless, because of the high ec- 
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centricity of the orbit and because of their adequate ob- 
servational coverage of the star near r.v. minimum, the 
final values of these three parameters are not expected 
to be significantly different. Assuming that the r.v. and 
the resolved companions are the same, we impose the 
additional constraints derived from our astrometric ob- 
servations, namely, the separation p and position angle 
between the binary components on 2003 May 13 (T = 
JD 2452772 = year 2003.36). Given the long (multi-year) 
orbital period and the small change (insignificant within 
the error bars) between our two relative astrometric ob- 
servations taken four months apart, we only use one of 
the astrometric measurements. Finally, we adopt a mass 
of 1. 05 M(7i for the GO V prim ary, based on an estimate 
from lDorren fc GuinanI l)1994|) . 

The equation that determines the binary orbit is Ke- 
pler's equation: 

S-esinS= ^(T-To), (3) 

where the eccentric anomaly E is related to the true 
anomaly 6 through 

tan-=(^^j tan-. (4) 
The remaining equations are: 

r = a(l — e cos E) (6) 
= r^{l-&iT?{e + uj)sm^i) (7) 

tan(^ — tan(0 -f cj) COS* 

tanS2 = - ; . (8) 

1 -I- tan (f) tan(6' + u) cos i 

Because the orbital period of HD 129333 is poorly con- 
strained by DM91, we choose to treat P as a free param- 
eter. Thus, the unknown parameters are eight: P, M2, 
i, a (the semi-major axis of the relative orbit), f2, E^ the 
radius vector r{T ~ Tq), and the true anomaly 9(T — Tq) 
of the companion in the relative orbit at time T — Tq. 
From the combined imaging and r.v. data we have im- 
posed seven constraints: p{T — Tq), (f>(T — Tq), e, uj. Mi, 
Ki , and T — Tq. Given that the number of unknowns is 
greater than the number of constraints, we cannot solve 
unambiguously for the parameters, let alone use a least- 
squares approach to determine their best-fitting values. 
However, by stepping through a grid of constant values 
for one of the parameters, we can determine the rest. 

We choose M2 as our step parameter for the grid, treat- 
ing it as a known parameter. In principle we can use 
Equations to express P in terms of M2, E, and the 
known variables, and then substitute this expression in 
Equation |3| which can be solved for E. However, be- 
cause of the complexity of the general functional form 
P{M2, E, p, (j), e,uj. Ml, Ki,T — Tq), and because Kepler's 
equation cannot be solved analytically, we employ a two- 
stage iterative approach. In the outer iteration, for a 
given value of M2 we converge upon a solution for P, 
and in the inner iteration we use the Newton-Raphson 
method to solve Kepler's equation for E. The iterative 
Newton-Raphso n method has be en described in detail 
elsewhere fe.g.. lPress et al.lll992|) . so we do not discuss 



it further. The convergence of the outer iteration loop 
however merits a brief description. 

We take an initial estimate Pq for the period and in- 
vert Equation |2] to obtain a numerical value for sin i as 
a function of P and M2 ■ From Equations we then 
express P in terms of the known parameters plus M2, E, 
and sini, and plug that expression in Equation |31 which 
is then easier to solve for E. Once E is obtained, we 
invert Kepler's equation to find a solution P = Pi for 
the orbital period that depends on the initial guess Pq. 
We repeat the above procedure by substituting Pi for 
Pq, and so forth until the values Pj converge. We stop 
when the value of P is constrained to better than 0.1%, 
which usually occurs after 3-4 iterations. Because of the 
monotonic dependence of the orbital elements i, a, and E 
on P, we can be certain that the thus-obtained solution 
for P is unique. 

Following the above procedure and adopting the DM91 
values for e, Tq, and Ki, we find that if the resolved 
companion is identical to the r.v. one, its mass is at least 
0.68Mq, with a corresponding period of 50.0 years, a = 
16.3 AU, and i = 85°. Values as small as M2 = 0.58Mq 
[P = 42.8 years, a = 14.4 AU, i = 84°) are possible if 
all parameters are set at their one-sigma deviations that 
minimize M2. 

However, a minimum mass of w O.58M0 for 
HD 129333B does not agree with the constraint from 
our IR spectroscopy, M2 < O.5OM0 (Sections lO 
obtained fr om comparison to t heoretical evolutionary 
tracks from iBaraffe et aJ] (|1998D . Moreover, a compan- 
ion with mass M2 > O.58M0 (spectral type K8 or ear- 
lier; iCoxl |200CI() wo uld be too br ight in absolute mag- 
nitude {Mk < 5.1: lBessel]lll991ii to be identified with 
HD 129333B (Mks = 6.30±0.12, from its apparent mag- 
nitude and from the Hipparcos distance to HD 129333). 
It is therefore likely that the r.v. and the spectroscopic 
companion are not identical. 

The inconsistency between the masses could be ex- 
plained b y noting that a recen t study of low-mass bi- 
naries by iHillenbrand fc White! l|2004l) has shown that 
most modern stellar evolutionary models tend to un- 
derestimate dynamical masses of main an d pre-main se- 
quence stars. According to the authors, the lBaraffe et alJ 
(1998) models underestimate the mass of a O.5M0 main 
sequence star by ~20%. This could reconcile the lim- 
its on the mass of HD 129333B obtained from near IR 
spectroscopy with those from the orbital solution. How- 
ever, the problem of the companion being sub-luminous 
remains. 

4.3.2. Comparison to Other Radial Velocity Data 

The DM91 set of r.v. data are the most accurate 
and deterministic for t he orb it of HD 1293 33. Other 
data ex ist from Wilson fc JoyI |T950,'l . Dorren fc GuinaiJ 
Ciaai, IMontes et a,].l ^2001al^. and from our own high- 
resolution spectroscopic observations^ (Figure El); how- 
ever, they do not improve the orbita l phase coverage 
greatly. Although iDorren fc GuinanI (data plotted as 
open squares) appear to have captured the binary near 
an r.v. maximum around 1993, their data are less restric- 

^ ''rad = -19.79 ± 0.37 km s-^ and -21.77 ± 0.62 km s^^ on 
2002 April 18, and 2003 February 10, respectively (Hillenbrand ot 
al., in prep.). 
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live because of the ir large uncertainties. In addition, the 
iDorren fc GuinanI data for 1990 systematically overesti- 
mate the r.v. of the primary with respect to DM91 data 
taken over the same period. W e therefore choose to dis- 
regard the IDorren fc GuinanI data set. The remaining 
data are very limited and we do not attempt to use them 
to re-fit the DM91 orbital solution. However, they are of 
some utility in constraining the orbital period. 

Because no other r.v. minimum is observed for 
HD 129333 between 1987.3 7 and 200 3.36, we conclude 
that P > 16 years. The Wilson fc Jovl data point (based 
on three measurements) is consistent with an r.v. min- 
imum, and is thus critical in constraining the orbital 
period. However, the auth ors do no t list an epoch for 
the observations. The Wils on~fc Jovl data were taken in 
the course of the Mt. Wilson stellar spectroscopic sur- 
vey, and are kept in the Ahmanson Foundation Star 
Plates Archive^. After consultation with the original 
plates, we retrieve the dates of the individual obser- 
vations: 1936 March 10 (year 1936.19), 1936 Jun 4 
(year 1936.42), and 1942 Jun 24 (year 1942.48). We 
adopt the mean date of these observations, the year 
1938.4l2:2: as the epoch for the lWilson fc Jovl measure- 
ment frad = —31. Oil. 3 km s~^, where the errors in the 
epoch correspond to the interval between their first and 
last observation. Because of its highly eccentric orbit, 
the star must have been within 1 year of r.v. minimum 
at this epoch. Given the r.v. minimum in 1987.37 and 
P > 16 years, we infer that the r.v. companion has a 
likely orbital period equal to the interval between the 
two observed minima, or to some integer fraction thereof: 
49.014:2 years, 2A.5t.2j years, or 16.31?^ years (aU con- 
sistent with the DM91 estimate). 

The 49-year orbital period agrees with the one ob- 
tained in Section 14.3.11 and supports the evidence that 
the r.v. and the resolved companion may be identi- 
cal. If the r.v. companion was on a 24.5-year period 
(a = 9.8 AU, from the DM91 orbital elements), it would 
have been > 0.33 — 0.44" from the primary during our 
imaging in 2003, with mass M2 > O.47M0. Such a 
companion should have been at least as bright as the 
resolved one {AKs = 3.0), although could have fallen 
just below our detection limits {AKg « 3.0 at 0.4") 
if it was at the lower limit of the allowed mass range. 
Given q = M2/M1 > 0.46 in this case, the star should 
be easily detectable as a double-lined spectroscopic bi- 
nary (SB2) through high-resolution spectroscopy in the 
near IR, where the contrast favors detecting SB2 systems 
with mass ratios as small as 0.2 llPrato e t al. 2002). A 
16.3-year period can most probably be excluded, since 
the 2003 data point does not indicate an approaching 
r.v. minimum (FigureUni) in late-2003-2004, as would be 
expected at this periodicity. 

Therefore, even after consideration of additional 
archival r.v. data, the question about the multiplicity 
of HD 129333 remains open. The system can be either 
a binary with a 14-16 AU semi-major axis (but a dis- 
crepancy in the inferred mass of the secondary), or a 
triple with a 10 AU inner (spectroscopic) and ~25 AU 
outer (resolved) components. Indeed, SIMBAD does list 
HD 129333 as a BY Dra variable, which may indicate 

Maintained at the Carnegie Observatories of Washington, 
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that the high level of chromospheric activity arises from 
close binarity, rather than extreme youth. However, the 
high photospheric Li I abundance of HD 129333 and its 
kinematic association with the Pleiades moving group 
(Section |4.2(l confirm its young age. Moreover, at 10 AU 
semi-major axis the inferred spectroscopic companion is 
too distant to be synchro nized with the rotation p eriod 
of the primary (2.7 days; |D orren fc GuinanI IT99|) . and 
to thus affect its chromospheric activity. An additional 
close-in component would be required, that should pro- 
duce a short-period SBl or SB2 spectroscopic signature, 
as in binary BY Dra systems. Such is not reported by 
DM91, however. 

In deciding whether a triple system with a 24.5-year 
period for the r.v. (inner) companion is a likely state for 
HD 129333, it is worth considering the dynamical sta- 
bility of such a system. We adopt masses of 1.05 Mq, 
0.5 Mq and 0.2 Mq for the primary, the inner, and the 
outer (resolved) companion, respectively, and apply a dy- 
namical stability criterion from the numerical analysis 
of Donnison & Mikulskis (1995). Assuming prograde or- 
bits, the|Donnison & Mikulskis condition for stability as 
applied to HD 129333 states that the distance of closest 
approach of the outer companion to the barycenter of the 
system should be >27 AU. Variations in the masses of 
the two companions within the determined limits do not 
change this requirement by more than 3-5 AU. At a pro- 
jected separation of 25.0 ± 1.5 AU from the primary, the 
resolved companion is fully consistent with this require- 
ment. Hence, the system can be a dynamically stable 
triple. 

4.4. HD 49197B: a Rare Young L Dwarf 

Our empirical knowledge of the photospheres of young 
(<1 Gyr) L dwarfs is currently very limited. The only 
confirmed such dwarfs are all companions to main se- 
quence stars: G 196-3B (L2, 20-300 Myr; Rebolo ct J] 
1998), Gl 417B5 (L4.5, 80-300 Myr; fKiritnatrick et all 
■2001.). the pair H D 130948B/C (L2, 300-600 Myr; 
Potter et all |2002|) . and now HD 49197B (L4, 260- 
790 Myr). It is useful to expand the sample of young 
L dwarfs in order to study gravity-sensitive features in 
their spectra, and to provide constraints for evolutionary 
models of ultra-cool dwarfs. 

Younger L dwarfs have bee n reported in several 
open clusters: a Ori (1-8 Myr- 'Zapatero Osorio et alJ 
[1999), the Trapezium ( ^1 Myr; [L ucas ct al. 2001), and 
Chameleon I (1-3 Mvr: iLopez Martf et al.ll2004t) . How- 
ever, these results have not been independently con- 
firmed. The youth of a Ori 47 (LI. 5), and hence its 
association with the cluster, has been r ecently brought 
into question bv iMcGovern et aD l)2004D . who find that 
the object shows strong K I absorption at J band, char- 
acteris tic of several Gyr old field dwarfs. iLucas et alJ 
l)2001l) determine Ml-lLs spectral types for their objects 
in the Trapezium, using water indices defined f or the 
i? fti 3 _ff-band spectra. They also use Burrows et alJ 
((1.997) sub-stellar evolutionary tracks to infer masses 
from IJH photometry. However, the deduced spectral 
types and the masses correlate very poorly — a result po- 

^ Gl 417B is itself considered to be resolved bv lBouv et alJ 12003) 
into two components with a 70 mas separation, equal to the diffrac- 
tion limit of their HST/WFPC2 observations. 



12 



Metchev & Hillenbrand 



tentially traceable to the anomalous continuum shapes of 
their spectra (their Figure 4), some of which appear to 
contain residual telluric or instrument-transmission fea- 
tures (as seen in their Figure 1) that the authors inter- 
pret as photospheric water absorption. Finally, in their 
analysis of photom etrically-identified brown dwarfs to- 
ward Chameleon I, iLopez Marti et all l)2004j) acknowl- 
edge that the classification of their early L dwarfs is 
uncertain, because their locus overlaps with that of ex- 
tincted M-type objects on their color-magnitude diagram 
(Figure 8 in that paper). 

Therefore, because of its association with a young star, 
HD 49197B is one of only five known young L dwarfs 
whose age can be determined with reasonable certainty. 
All five span a narrow range in spectral type: L2-L5. A 
program of uniform spectroscopic observations of these 
young L dwarf companions, undertaken in a manner sim- 
ila r to the NIRSPEC b rown dwarf spectroscopic survey 
of IMcLean et al.l lj2003(l . p romises to establish grav ity- 
sensitive standards (as in IKleinmann fc Halll [T98^ . for 
F8-M7 stars) to use in determining the ages of isolated 
L dwarfs. 

4.5. Sub-Stellar Companions to Main-Sequence Stars 

Until recently, only a handful of brown dwarf com- 
panions to nearby A~M stars were known from di- 
rect imagin g, all at angular separations >4" (see com- 
pilation in iReid et J] |200H) . With AO technology 
still in its early developing stages, ground-based direct 
imaging observations of main-sequence stars were sen- 
sitive mostly to massive, widely-separated sub-stellar 
companions. From the observed dearth of brown 
dwarf companions to main sequence stars at sepa- 
rations comparable to those in main sequence bina- 
ries, it was inferred that the radial-velocity "brown 
dwarf de sert" (for sep arations <3 AU: Campbell et al.l 
liSSS: Ma rcv fc BeTHta flQSQ: Marcv & Butler 2000) ex- 
tended to at least 120 AT J ^pp penhcim er et al...200Ll. 
or 12 00 AU l)McCarthvll200in McCarth v fc ZuckermanI 
|2004|) . with estimates for the brown dwarf companion 
frequency around 1% wit hin th is sepa ratiori range. From 
2MASS data however, ICizis et alJ l)200lD found that 
the brown dwarf companion fraction was much higher 
(~18%) at separations >1000 AU from F-MO dwarfs, 
and consistent with that of stellar companions to G stars 
((puaucnnov & Mavor 1991). Though the Gizis ct al] re- 
sult is based on only 3 bound companions out of 57 then 
known field L and T dwarfs (a fourth bound companion, 
G J 1048B, has now been confirmed in the same sample 
bv lSeifahrt. Neuhauser. fc MugraueJ2004(l . they exclude 
a brown dwarf companion fraction of 1.5% at the 99.5% 
confidence level. Such an abrupt change in the frequency 
of bound brown dwarfs at 1000 AU from main-sequence 
stars is not predicted by any of the current brown dwarf 
formation scenarios. More likely would be a continu- 
ously varying sub-stellar companion fraction from inside 
the r.v. brown dwarf desert at <3 AU out to distances 
>1000 AU. 

Recent results from more sensitive space- and ground- 
based surveys point to a somewhat higher frequency 
of sub-stellar companions. In a survey of 45 young 
stars within ~50 pc, the NICMOS Environments of 
Nearby Stars team has reported t he discovery of 2 
confirmed brown dwarfs, TWA 5B l)Webb et alJ 119991: 



i Lowrance et all Il999|) and HR 7329B ijLowrance et alJ 
120001) ■ and a probable third one: the binary com- 
panion Gl 577B/C, whose compone nts likely span the 
stellar/substellar boundary (Lowran ce et al.l 12003(1 . A 
similar survey of twenty-four 5-15 Myr old stars in 
the more distant (wl50 pc) Scorpius-Centaurus asso- 
ciation d oes not detect any pl ausible sub-stellar com- 
panions l|Brandner et alJ l2000|) . Even so, the frac- 
tion of stars with sub-stellar companions detected with 
NICMOS (2-3 out of 69) is markedly higher than the 
one reported from the two initial large-scale ground- 
based surveys (2 out of »390: [ Qpnenheimer et alJ 
l200lt iMcCarthv fc ZuckermanI 120041) . and is inconsis- 
te nt with the in completenes s-corr ected estimate (<2%) 
of iMcCarthv fc ZuckermanI l|2004|) . Furthermore, as a 
result of improvements in existing AO technology, and 
the equipping of several large telescopes with newly- 
designed high-order AO systems, recent ground-based 
direct imaging efforts have been more successful in de- 
tecting close-in bro wn dwa rf companions to sun-like 
primaries: Gl 86B l|Els et a l. 2 001), HP 130948B/C 
i|Potter et al1l20^ . HR 7672B l(Liu et al J 120021 and 
HD 49197B (this paper). All of these are at angular sep- 
arations <3", and at projected distances <50 AU from 
their primaries, and hence inaccessible for imaging by 
McCarthy & Zuckerman (2004), whose survey targeted 
the 75-1200 AU separation range. Finally, with the nat- 
ural guide star limit of AO systems being pushed to 
ever fainter magnitud es using curvature sensors (down 
to ~16 mag at 0.8um: 'Sieerl er et al]l2003D . a number of 
very low-mass (VLM) binaries has become known, with 
separations as small as 1 AU. The components in these 
VLM binaries often straddle the stella r-substellar b ound- 
ary (for a compilation, see Table 4 in lClose et al.l l2003'). 

The emergent picture from these recent discoveries 
is that of potential deficiency of brown dwarfs at 10- 
1000 AU separations from main sequence stars, though 
not as pronounced as in the r.v. brown dwarf desert 
(frequency <0.5%: IMarcv fc ButleJ l2000|) . Based on 
one detection (o f a binary brown dwarf companion) 
among 31 stars, IPotter et al . (2002) set a lower limit 
of 3.2±3.2% for the frequency of brown dwarfs at 10- 
100 AU from main sequence stars. At separations 
>50 AU, fr om the NICMOS discoveries ijLowrance et alJ 
'200(f) and from their newly-reported brown dwarf 
com panio n to the star GSC 08 047-00232 in Horologium, 
iNeuhauser fc Guentheil l)2004D report that brown dwarfs 
are found around 6ib4% of sta rs. The outer scale for the 
INeuhauser fc Guentheil l|2004f) estimate is not specified, 
but is probably Hmited to 1000-2000 AU by the FOV 
of high angular resolution IR detectors (up to 40"x40"; 
e.g., NICMOS, or ones used behind AO), and by the dis- 
tances out to which young stars are probed for sub-stellar 
companions (out to 100-200 pc). B y combining these es- 
timates with thel Gizis et al.l l(200lD estimate of 18±14% 
at separations >1000 AU, we can conclude that, despite 
the small number statistics involved, there possibly exists 
a continuum in the frequency distribution of brown dwarf 
companions at separations ranging from within the r.v. 
brown dwarf desert (<3 AU) out to >1000 AU (where 
brown dwarf companions are as common as stellar ones). 
The observed decline in the rate of occurrence of directly 
imaged brown dwarf companions at small separations is 
likely at least partially an effect of the limited sensitivity 
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of imaging surveys to close-in low-mass brown dwarfs. 
New, sensitive surveys for sub-stellar companions, such 
as the Palomar AO Survey of Young Stars are poised to 
explore this regime in the next few years. 



5. CONCLUSION 

We have presented the observing strategy and first 
results from the Palomar Adaptive Optics Survey of 
Young Stars, aimed at detecting sub-stellar companions 
to <400 Myr solar analogs within 160 pc of the Sun. We 
have discovered low-mass (0.04-0.5 Mq) companions to 4 
young nearby stars. The L4±l brown dwarf HD 49197B 
and the M2±l V star HD 129333B have been confirmed 
as companions to their corresponding primaries through 
follow-up astrometry and spectroscopy. Physical associ- 
ation in the V522 Per and RX J0329.1-t-0118 systems, 
containing respectively M4±l and M3±2 secondaries, 
has been established with >99.95% confidence from spec- 
troscopy and from the expected space density of objects 
of similar spectral type. 

The astrometry for the resolved stellar companion to 
HD 129333 is found to be consistent with archival r.v. 
data for this single-lined spectroscopic binary, indicat- 
ing that the resolved and the r.v. companions may be 
identical. Given the inferred mass constraints on the 
secondary however, the companion is then under lumi- 
nous by at least 1 mag at Ks- A solution in which 
the star is a triple is also likely. It does not suffer from 
similar inconsistencies, and could be dynamically stable. 
Because the expected mass ratio between the inner two 
companions of the triple is >0.46, they should be resolved 
as a double-lined spectroscopic binary from high resolu- 
tion infrared spectroscopy. In either case HD 129333 is 
confirmed to be a multiple star, and hence not a true 
analog of the y oung su n, as previously co nsidered (e.g., 
IDorren fc Cu Tnan 1994 IStrassmeier fc RTce 1 998^ . 

The newly-discovered brown dwarf HD 49197B is one 
of very few confirmed young (<1 Gyr) L dwarfs. It is also 
a member of a small number of brown dwarf companions 
imaged at projected separations of <50 AU from their 
host stars, i.e., at distances comparable to the giant- 
planet zone in the Solar System. The number of such 
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companions, albeit small, has been growing steadily in 
recent years with the results of more sensitive imaging 
surveys coming on-line. Longer duration radial velocity 
surveys and improvements in AO techniques are expected 
to further push the detection limits of each method to the 
point where their sensitivities overlap. Although the true 
extent and depth of the so-called "brown dwarf desert" 
will not be revealed until that time, increased sensitivity 
to sub-stellar companions at small separation has already 
resulted in upward revisions of their estimated frequency. 
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TABLE 1 
Observations. 



Target 


First Epoch 


Mode 


Second Epoch 


Mode 




Telescope 




Telescope 




HD 49197 


2002 Feb 28 


JHKg coronagraphic imaging 


2003 Nov 9-10 


JHKs coronagraphic imaging. 




Palomar 




Keck II 


K spectroscopy 


HD 129333 


2003 Jan 12 


JHKs non-coronagraphic imaging, 


2003 May 13 


Br7 non-coronagraphic imaging 




Palomar 


JK spectroscopy 


Palomar 




V522 Per 


2003 Sep 20 


JHKs non-coronagraphic imaging 


2003 Nov 10 


KsL' non-coronagraphic imaging. 




Palomar 




Keck II 


JK spectroscopy 


RX J0329.1+0118 


2003 Sep 21 


JHKs coronagraphic imaging 


2003 Nov 10 


K spectroscopy 




Palomar 




Keck II 





TABLE 2 

Properties of the observed stars. 



Object 


P.M. (mas 

lla COS 6 


/year) 


Parallax 
(mas) 


Spectral type 


Ks'' 
(mag) 


Age'' 
(Myr) 


Notes 


HD 49197 

HD 129333 

V522 Per 

RX J0329.H-0118 


-35.12 ± 1.05 
-138.61 ±0.61 
17.6 ±3.0 
5.4 ± 1.1 


-48.59 ± 0.63 
-11.92 ± 0.68 
-26.9 ±2.7 
-5.8 ± 1.1 


22.41 ±0.87 
29.46 ±0.61 
5.46 ±0.20 
< 10 


F5 V 
GO V 

GO IV 


6.067 ±0.024 
5.914 ±0.021 
9.352 ± 0.024 
9.916 ±0.019 


260 - 790 
10 - 100 
90 ± 10 
120 


1,3 
1,4 
2,5 
2,6,7 



''From the 2MASS Point Source Catalog 

''See Sectional 

Notes: 1. Proper motion and parallax from Hipparcos IPerrvman et alj|1997). 2. Proper motion from Tycho 2 iH0g et al ■■20001) ■ 
3. Spectral type from Hiyy arcos ^Ferry man ct al. 1.9^2 t) • 4. Spectral type from Buscombe Fosteil il997D . 5. Assumed to be at the 
mean Hippa rcos distance (v an Leeuwerill999i) of the a Per cluster. 6. A distance of at least 100 pc can be inferred for this young 
(~100 MyrjFVink ct al. 19 ^3) star from its small proper motion and its location toward the Taurus star forming region. 7. Spectral 
type from lBuscombel fig 



TABLE 3 

IR MAGNITUDES AND COLORS OF THE COMPANIONS. 



Object A J AH AKs J -H H - Ks Ks Ks - L' 



HD 49197B 


9.6 


± 


1.2 


8.52 


± 


0.12 


8.22 


± 


0.11 


1.2 


± 


1.2 


0.33 


± 


0.20 


14.29 ±0.14 


HD 49197"C" 


6.86 


± 


0.10 


6.82 


± 


0.09 


6.68 


± 


0.10 


0.27 


± 


0.14 


0.17 


± 


0.14 


12.75 ±0.10 


HD 129333B 


3.38 


± 


0.10 


3.13 


± 


0.09 


3.04 


± 


0.08 


0.55 


± 


0.14 


0.19 


± 


0.12 


8.95 ±0.08 


V522 PerB 


5.69 


± 


0.09 


5.44 


± 


0.09 


5.16 


± 


0.09 


0.63 


± 


0.13 


0.38 


± 


0.13 


14.51 ±0.09 0.15 ±0.18 


RX J0329.1+0118B 


4.22 


± 


0.12 


3.86 


± 


0.08 


3.65 


± 


0.08 


0.59 


± 


0.15 


0.31 


± 


0.13 


12.85 ±0.09 



TABLE 4 

ASTROMETRY OF THE COMPANIONS. 



Object Epoch 1 

offset (arcsec) P.A. (degrees) 



Epoch 2 Epoch 1 (if non-common p.m.) 

offset (arcsec) P.A. (degrees) offset (arcsec) P.A. (degrees) 



HD 49197B 

HD 49197"C" 

HD 129333B 

V552 PerB 

RX J0329.1+0118B 



0.9499±0.0054 

6.971±0.030 
0.7343±0.0032 
2.0970±0.0090 
3.75±0.05 



78.25±0.40 
346.13±0.34 
173.19±0.35 
194.02±0.34 
303±5 



0.9475±0.0022 

7.016±0.008 
0.7363±0.0032 
2.0937±0.0032 
3.781±0.016 



77.60±0.25 
346.50±0.10 
173.37±0.35 
193.91±0.11 
303.85±0.34 



0.9029±0.0037 

6.950±0.009 
0.7221±0.0033 
2.0991±0.0032 
3.714±0.070 



81.87±0.31 
346.10±0.10 
180.68±0.37 
193.93±0.11 

303.7±1.1 
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TABLE 5 

Spectroscopic measurements for the companions. 



Star 


EW{Na I)'' 

(A) 


EW(Ca 1)=^ 

(A) 


EW(CO)=' 

(A) 


EW(K I)'' 

(A) 


H2O index'" 


(K) 


HD 49197B 


-0.3 ±0.6 


0.2 ±0.4 


0.82 ± 0.05<* 






<3000 


HD 129333B 


5.17 ±0.59 


4.79 ± 0.64 


7.26 ±0.57 


2.18 ± 0.11 




3660 


V522 PerB 


5.87 ±0.58 


3.11 ±0.50 


5.67 ±0.38 


1.07 ±0.16 


0.91 ±0.01 


3200 


RX J0329.1+0118B 


6.76 ±0.60 


3.81 ±0.39 


8.63 ±0.17 






3300 



^At K band as defined bv lAli et all 1199511 . 

''At J band as defined bv lGorlova et alj 1200311 . 

'^Calculated from empirical relations relating T^ff to EW (Ca I) (Table 4 in lAli et alj|199'5ll . The sign of the 
linear coefficient in the Ca I "cool" relation of lAli et alj has be en change d from ' (as erroneously listed in 
their paper) to ' + ' to match the slope of their empirical data. lAli et all quote an error of ±300 K for this 
index. 

''The CO index for HD 49197B is not in A, but as defined bv lMcLean et alj 1200311 . 



TABLE 6 

Estimated properties of the companions. 



Objects 


Spectral type 


Mass 


Projected Separation 
(AU) 


Probability of 
chance alignment 


HD 49197B 

HD 129333B 

V522 PerB 

RX J0329.1+0118B 


L4 ± 1^ 
M2 ± 1'' 
M4 ± 1*" 
M3 ± 2*" 


n 060+°"i2 

20+°'3° 
0.125 ± 0.025 
n 20+"-^° 


43 
25 
400 

380 


3 X lO"'' 
3 X lO"** 
2 X 10-3 
9 X 10-* 



^Inferred from the iC-band spectrum, and from the absolute magnitude of the object fSection l3.3.1i . 
Based on the estimate of the Ca I-derived effective temperature (Table and on the J-band K I 
and H2O abso rption (if available). A MK spectral type vs. Tcff classification for dwarfs is adopted from 
IBessell 1199111 . 
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Fig. 1. — A'5-band (2.15/im) first (a, b) and second (c) epoch images of HD 49197. Panel (a) shows both candidate companions to 
HD 49197; panels (b) and (c) arc zoomed in to point out only the close-in bona-fide companion (HD 49197B). The first-epoch image is the 
result of 24 median-combined 60 sec exposures with Palomar/PHARO, whereas the second-epoch image was formed by median-combining 
six 60 sec exposures with Keck/NIRC2. A 1.0"-diameter coronagraph occults the primaxy in both cases; in the Keck image the coronagraph 
shows a residual ft;0.16% transmission. HD 49197B was initially unnoticed in the first-epoch image, where its detection was hindered by 
the presence of equally prominent AO speckles. 



HD 138333 VSE2Per RX J0320.1+oi.ia;: 



■fi 




Fig. 2. — Palomar images of the stellar companions. The HD 129333 image is taken through a narrow-band (1%) Br7 (2.166 fj,m) filter, 
while the V522 Per and RX J0329.H-0118 images are taJcen at Kg. Five dithered 1.4 sec exposures were aligned and median-combined to 
obtain each of the displayed images. 
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Fig. 3.— K-band spectra of the low-mass companions from Palomar (HD 129333B; R « 1000), and Keck (RX J0329.1+0118B, V522 PerB, 
and HD 49197B; R R! 2700). All spectra have been normalized to unity at 2.20/um and offset by 0.5 in the vertical axis. 
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Fig. 4.— J-band spectra of two of the low-mass companions from Palomar (HD 129333B; R 1200), and Keck (V522 PerB; R « 2400 
in the fifth order, and R Ri 2900 in the sixth). All spectra have been normalized to unity at 1.25;um and offset by 0.5 in the vertical axis. 
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Fig. 5. — Near-IR color-color diagram of the detected companions. The solid line represents the B2V— L8V main sequence, with data 
compiled from Cox (2000, B2V -M6y) and Kirkpatrick ct al. (2000, M8V-L8V). The d ashed line shows the G 0-M7 giant branch ICmJ 
|200y). The fCox. and the K irkpatrick et aH colors are converted from the Johnson-Glass fBessell & Bret1?'1988} and 2MASS ifnutrTe^an 
12003) systems, respectively, to the CIT system using relations from 'Carpenter (2001). HD 49197"C" is too blue to be a bona-fide low-mass 
companion to HD 49197. The near IR colors of the remainder of the companions agree well with their inferred spectral types (Table IB! . 
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Fig. 6. — Proper motion diagram for tlic two companions to HD 49197. The ofTscts from the primary at cacfi observational epocfi arc 
plotted as solid points with errorbars. The inferred offsets at the first epoch (assuming non-common proper motion) are shown just with 
errorbars. The dashed lines reflect the proper and parallactic motions of the primary between the two epochs. HD 49197B is a common 
proper motion companion within the 1 sigma errors, while HD 49197"C" is more consistent with being a ba<;kground object. 
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Fig. 7. — Same as Figure|n]for the companion to HD 129333. Within the 1 sigma errors, the companion shares the proper motion of the 
primary. 
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Fig. 8.— Same as FigurelHlfor the companions to V522 Per (left) and RX J0329.1-f 0118 (right). The Nov 21, 1995 data point for the 
companion to RX J0329. 1-1-0118 is from Sterzik et al. 11997) and is set to be at the mean epoch of their observations (Nov 19-23, 1995). 
The proper motion of the primaries between the observations epochs are too small to decide the physical association of th e com panions 
within the astrometric errors. The probability of association in each system is estimated from near IR spectroscopy ("Section 14. II . 



Palomar AO Survey of Young Stars 



25 



10 



a 
o 
o 
m 
O 
u 

o 

CD 

Oh 



o 5 







"1 1 r 



1 1 1 r 



1 1 r 



HD 49197« 

15 fi/ 



RX J0329.1+0118B 

/ 

V522 PerB 



-r / 
/ 



/ 

/ 



/ 

HD /29333B 



^ 



HD 49197"C" 



J I I L 



J I I L 



J I L 







5 10 
from photometry 



Fig. 9. — Comparison of the photometrically derived absolute Xg-band magnitudes of the companions (assuming the hehocentric 
distances of the corresponding primaries) to the spectroscopically inferred ones. The dashed line has a slope of unity. The location of 
HD 491 97"C " (of which no spectra were taken) along the vertical axis is based on a spectral type (F— G) inferred from its near IR colors 
fSection l3.1l . HD 49197"C" is intrinsically too bright to be at the same heliocentric distance as HD 49197, whereas the remainder of the 
companions are consistent with being at the distances of their respective primaries. 



26 



Metchev & Hillenbrand 




1940 1960 1980 2000 

Year 



Fig. 10. — Radial velocity data for HD 129333 from the literature. Data from DM91 arc plotted as filled circles, from'Dorren & GuinaiJ 
as open squares, from Monies ct al. 1 2001a) as filled sq uares, an d from Wilson & Joy 1 1950!) as an open circle. Our own unpublished 
r.v. data are shown as solid triangles. The exact epoch of the lWilson Jovi observation is unknown. Given the DM91 orbital solution and 
the current phase coverage, periods <16 years can be excluded. 



